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Abstract 

While over the last century or more considerable effort has been put into 
the problem of finding approximate solutions for wave equations in general, and 
quantum mechanical problems in particular, it appears that as yet relatively 
little work seems to have been put into the complementary problem of estab- 
lishing rigourous bounds on the exact solutions. We have in mind either bounds 
on parametric amplification and the related quantum phenomenon of particle 
production (as encoded in the Bogoliubov coefficients), or bounds on trans- 
mission and reflection coefficients. Modifying and streamlining an approach 
developed by one of the present authors [Phys. Rev. A 59 (1999) 427-438], 
we investigate this question by developing a formal but exact solution for the 
appropriate second-order linear ODE in terms of a time-ordered exponential of 
2x2 matrices, then relating the Bogoliubov coefficients to certain invariants 
of this matrix. By bounding the matrix in an appropriate manner, we can 
thereby bound the Bogoliubov coefficients. 

Keywords: Bogoliubov coefficients, Transmission coefficient, 

Reflection coefficient, rigorous bounds. 
larXiv:0801.06T0V l [quant-ph] 
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1 Introduction 



There are numerous physical situations in which it is both extremely interesting and 
important to study the second-order ODE p] 

4>(t) + u 2 (t)<p(t) = o, (l) 

or its equivalent in the space domain [I] 

4>"{x) + k 2 (x)(f)(x) = 0. (2) 

Viewed in terms of the time domain, equation (0Q) can be viewed as an example 
of parametrically excited oscillation; it arises for instance when a wave propagates 
through a medium whose refractive index is externally controlled to be a function 
of time (though remaining spatially invariant) Q In contrast, the spatial version of 
this equation as presented in ([2]) arises classically in situations where the refractive 
index is spatially dependent (so called "index gradient" situations), or in a quantum 
physics context when considering the Schrodinger equation for a time-independent 
potential: 

- lf-<l>"(x) + V(x)<l>(x) = E<j>(x), (3) 
2m 

as long as one makes the translation 

;2 , 2m[E-V{x)} fA . 
k 2 (x) ~ h2 WJ . (4) 

However they arise, equations ([1]) and ([2]) are central to the study of both quantum 
physics and wave phenomena generally. 

Because of this central importance, over the last century or more a vast body 
of work has gone into the question of finding approximate solutions to equations 
(PQ) and (j2J), most typically based on JWKB techniques and their variants (phase 
integral techniques, etc.) [5]. In contrast very little work seems to have gone into the 
physically important question of finding explicit bounds on the relevant Bogoliubov 
coefficients and/or reflection and transmission coefficients [TJ. 

In the current article we shall modify and streamline the analysis of [1] ; presenting 
an alternative proof that is considerably more direct and focussed than that in [1] . 



^or instance, situations of this type have been used to model sonolumincscence [2J, and more 
recently both quasiparticle production in analogue spacetimes [3] and analogue signature change 
events [4]. In all these situations it is extremely useful to have rigorous and largely model- 
independent bounds on the amount of particle production that might reasonably be expected. 
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To keep the discussion simple and straightforward we shall assume that oj(t) — > ujq 
(equivalently k(x) — > k ) outside some region of compact support [tj,t/] (equiva- 
lently [xi,Xf]). That is, concentrating on the time-domain formulation of equation 
(JT]), the quantity co 2 (t) — cu 2 is a function of compact support H Because of this com- 
pact support property we know that everywhere outside the region the exact 
solution of the wave equation (CQ) is given by linear combinations of exp(±iu;oi), 
and that the central question to be investigated is the manner in which exact so- 
lutions on the initial domain (—00, ti) "connect" with exact solutions on the final 
domain (tj,+oo). Describing and characterizing this "connection" is exactly what 
the Bogoliubov coefficients are designed to do. 



2 Time-ordered exponentials 

We are interested in solving, exactly but possibly formally, the second-order PDE 

]>(t) +uj 2 (t)(f)(t) = 0. (5) 
One way of proceeding is as follows: Define a momentum 

7T = 0, (6) 

and then rewrite the second-order ODE as a system of first-order ODEs 

= tt; (7) 

vr = -uo 2 (t) 0; (8) 

or in matrix notation (where we have carefully arranged all matrix elements and 
vector components to carry the same engineering dimensions) 



d 







dt 














_ 7r/^0 _ 



(9) 



u 
-u 2 /uj 

This matrix ODE always has a formal solution in terms of the so-called "time ordered 
exponential" 



J t 



2? < exp 



to L 



0_ ^0 
-u 2 {t)/u)Q 



dt 



tt/uj 



(10) 



J to 



2 This "compact support" condition is not strictly necessary, and at the cost of a little more 
analysis one can straightforwardly extend the comments below to a situation where there is a finite 
limit u>(t) — > ujoo as t — > ±00 PQ. At the cost of somewhat more tedious additional work, there 
are also useful things that can be said of the situation where uj(t) — > w±oo, with (x>_oo 7^ lo +00 , as 
£-►±00 p]. 
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The meaning of the time-ordered exponential is somewhat tricky, but ultimately is 
just a 2 x 2 matrix specialization of the operator-valued version of the "time ordered 
exponential" familiar from developing quantum field theoretic perturbation theory 
in the so-called "interaction picture" [6]. Specifically, let us partition the interval 
(t , t) as follows: 



t < ti < t 2 < t$... < t„_3 < t n _2 < t n -i < t n — t, 
and define the "mesh" as 

M = max {ti — ti-i}. 

ie(l,n) 

Then define the time-ordered exponential as 



12) 



T(t) 



ST { exp 



to 

n-1 



u 
-u 2 (t)/u 



dt 



M->0, (r. 



lim [ exp 

0. (n-»oo) A A 
i=0 







-u 2 (t n -i)/uj 



(t, 



■"n—i—X, 



(13) 



Note that in this matrix product "late times" are always ordered to the left, and 
"early times" to the right. By working with this time-ordered matrix we will be 
able to extract all the interesting physics. (If we work in the space domain then 
the equivalent matrix T is "path-ordered", and is closely related to the so-called 
"transfer matrix".) 

• Since all of the "complicated" physics takes place for t £ (ti, tf), it is also useful 
to define 



T = & 



( exp (/,' 



u 
-u) 2 (t)/u)Q 



dt 



a b 
c d 



(14) 



We are guaranteed that det[T] = 1, that is ad — be — 1. This follows from the 
fact that det[T] = exp{tr(ln[T])}, and the explicit formula for T above. 

Another particularly nice feature is that with the current definitions the trans- 
fer matrix T is manifestly real. This is relatively rare when setting up scattering 
or particle production problems, so we shall make the most of it. 
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3 Bogoliubov coefficients 

Let is now calculate the Bogoliubov coefficients. Before ti, and after tf, the wave- 
function is just linear combinations of exp(±iu t). We can prepare things so that 
before ti the wavefunction is pure exp(+iuot), 

ip(t < ti) = exp(+iu t); (15) 

in which case after tf the wavefunction will be a linear combination 

^(t > tf) = aexp(+iu t) + /3exp(— iu t), (16) 

where the Bogoliubov coefficients a and (3 are to be calculated. That is, we have 

(17) 








exp(-Hco>o ti) 


_ ?tM) _ 




i exp(+iu> 1^ 



and 



But we also have 













aexp(+iuj tf) + /3exp(— iu tf) 
i {a exp(+iuotf) — (3 exp(—iuj tf)} 



T 



TX/UJQ 



(18) 



(19) 



implying 



a exp(+ia;o tf) + (3 exp(— iu tf) 
i {a exp(+iu tf) — (3 exp(— iu tf)} 



a exp(+io;o U) +bi exp(+iu ti) 
cexp(+iu ti) + diexp(+iu ti) 



Solving these simultaneous linear equations we find 

1 



a 



ft 



a + d + i(b - c)] exp(-iu [tf — ti]), 
a - d + i(b + c)] exp(-iuJo [tf + U]), 



(20) 



(21) 



(22) 



so that the Bogoliubov coefficients are simple linear combinations of elements of the 
matrix T. Then (remember the matrix T is real) 



|a| 2 = i {(a + rf) 2 + (6-c) 2 } 



(23) 
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\P\ 2 = \{(a-df + (b + cf}, (24) 

and so 

| ,2 | fl ,2 (g + rf) 2 + (6 - C) 2 - (O - rf) 2 - (b + c) 2 

2ad - 2bc + 2ad - 26c , , „ . . 

= —ad — 00= 1, (26) 

thus verifying that, (thanks to the unit determinant condition), the Bogoliubov coef- 
ficients are properly normalized. Particle production is governed by the (3 coefficient 
in the combination 

|/?| 2 = i{(a-d) 2 + (6 + c) 2 }, (27) 

= - {a 2 + d 2 - 2ad + b 2 + c 2 + 2bc] , (28) 

= I{a 2 + rf 2 + 6 2 + c 2 -2}, (29) 

= itr{TT T -/}. (30) 

Note that the transpose T T is now time- anti- ordered. 
Similarly 

\a\ 2 = i{(a + d) 2 + (6-c) 2 }, (31) 

= ^{a 2 + d 2 + 2ad + b 2 + c 2 -2bc} , (32) 

= ^{a 2 + d 2 + b 2 + c 2 + 2} , (33) 

= ^tr{TT T + /}. (34) 

In summary, we can always formally solve the relevant ODE, either equation (CQ) or 
its equivalent equation (j2D, in terms of the time-ordered exponential, and we can 
always formally extract the Bogoliubov coefficients in terms of traces of the form 
tr{TT T }. We shall now use these formal results to derive rigorous bounds on the 
Bogoliubov coefficients. 
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4 Elementary bound: 

Consider the quantity 



X(t) = T(t) T(t) T = ST 



x 



f xp (Z 



& { exp 



oo 
-uj 2 {t)/u) 



dt 



u 
-u 2 (i)/uj 



di 



(35) 



dX 
~dt 





-cj 2 (t)/cu 



This object satisfies the differential equation 

X(t) + X(t) 
with the boundary condition 

Now note 



-cu 2 (t)/cuo 
u; 



X(t t ) = /. 
tr(X) = tr{TT r } = a 2 + b 2 + c 2 + d 2 . 



(36) 

(37) 
(38) 



Furthermore 
dX 

~dT 



a 2 + b 2 ac + bd 
ac + bd c 2 + d 2 

a 2 + b 2 ac + bd —uj 2 /ujq 
ac + bd c 2 + d 2 uj 

2oo (ac + bd) u (c 2 + d 2 ) - (uo 2 /uj )(a 2 + b 2 ) 

u (c 2 + d 2 ) - (u 2 /u )(a 2 + b 2 ) (-2LU 2 /u )(ac + bd) 



u 
-uj 2 /uj 



+ 



(39) 



and so we see 
tr 



uj 
uj 2 /u)q 



Therefore 



X + X 



dtx[X] 
dt 



—uj 2 /uj 
u 



2(ac + bd) 



2(ac + bd) 



u 



^0 



UJ 

UJ 



(40) 



(41) 



Using this key result, and some very simple analysis, we shall now derive our first 
elementary bound on the Bogoliubov coefficients. 
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For any 2 real numbers, using (x + y) 2 > and (x — y) 2 > 0, we have 



x 2 + y 2 > 2|xy| 



In particular, for any 4 real numbers this implies 

^ 2 + b 2 + c 2 + d 2 > 2y/(a 2 + b 2 )(c 2 + d 2 ). 



a 



(42) 



(43) 



But we also have 



\ac + bd\ 2 + \ad - bc\ 2 = a 2 c 2 + 2abcd + b 2 d 2 + a 2 d 2 - 2abcd + b 2 c 2 (44) 

= (a 2 + b 2 )(c 2 + d 2 ), (45) 



thus, for any 4 real numbers 



+ b 2 + c 2 + d 2 > 2 v /\ac + bd\ 2 +\ad-bc\ 2 . 



(46) 



For the particular case we are interested in we additionally have the unit de- 
terminant condition ad — be = 1, so the above implies 



whence 



a 2 + b 2 + c 2 + d 2 > 2^/\ac + bd\ 2 + 1, 
2\ac + bd\ < v/(a 2 + b 2 + c 2 + d 2 ) 2 - 4. 



Then 

whence 

dtrjXj 
dt 

whence 



dtr[X] 
dt 



2(ac + M) 



a; 

u; 

w 



< 2\ac + bd\ 



w 

^0 



< y/(a? + b 2 + c 2 + d 2 ) 2 - 4 



w 

w 

W 



This implies 



1 dtrjX] 
^tr[X] 2 - 4 dt 

dcosh -1 tr[X/2] 



< 



tr[X] 2 - 4 



W 



a; 



w 



w 
w 



df 



< 



w 

w 

w 



(47) 
(48) 

(49) 

(50) 
(51) 
(52) 
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whence 



We now have 





uo 2 


dtj 


{f 


LUq 




iO 





tr[X] < 2 cosh 



\P\ 2 = ]{tr{TT T }-2} = 1 -{tr{X}-2} 



so that 



\P\ 2 < 2^ 



u 

UJ 



dt 







u 2 


dtj 


{i 




u 






^0 





So finally 



and consequently 



\(5\ 2 < sinh 2 



la; I < cosh < - 









u 2 


dtj 




r 


L)Q - 


















u 2 


dtj 




f" 


LU - 






hi 




UJ 





(53) 
(54) 

(55) 
(56) 

(57) 
(58) 



These bounds are quite remarkable in their generality. A version of this result was 
derived in [1] but the present derivation is largely independent and has the virtue 
of being utterly elementary — in particular, the use of complex numbers has been 
minimized, and we have completely eliminated the use of the "auxiliary functions" 
and "gauge conditions" that were needed for the derivation in [TJ . 

If one translates this to the space domain, then the equivalent barrier penetration 
coefficient is T transmission <-» l/|a| 2 , and the equivalent reflection coefficient is R <-» 
|/3 2 |/|a| 2 . Making the appropriate translations 



1] 



transmission 



> sech 





k 2 (x) 


dx| 




k 







and 



R < tanh 2 



k 2 {x) 



dx 



(59) 



(60) 



(For completeness we mention that reference [T] provides a number of consistency 
checks on these bounds by comparing them with known exact results [7].) 
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5 Lower bound on \(3\ 2 

To obtain a lower bound on the \(3\ Bogoliubov coefficient, consider any real valued 
parameter e. Then since the matrix T is itself real, 

tr {(T - eT T ) T (T - eT T )} > 0, (61) 

so that 

(1 + e 2 ) tr(T T T ) - 2e tr(T 2 ) > 0, (62) 

whence 

tr(T r T) > tr(T 2 ), (63) 

1 + e A 

This bound is extremized for e = ±1, whence 

tr(T T T) > |tr(T 2 )| , (64) 

and so 

|/3| 2 >i{|tr(T 2 )|-2}. (65) 

This is certainly a bound, but it is not as useful as one might hope. It is useful only 
if tr[T 2 ] > 2. But 

tr[T 2 ] = a 2 + d 2 + 2bc = a 2 + d 2 + 2(ad - 1) = (a + d) 2 - 2 = (tr[T]) 2 - 2. (66) 

So using the unit determinant condition, tr[T 2 ] > 2 can be seen to require |a+d| > 2, 
that is, tr[T] > 2. But when does this happen? For the real matrix 



a b 
c d 



(67) 



with unit determinant the eigenvalues are 



x=a _ ±1±V(a + dr -4 (6g) 

The condition a + d > 2 is thus equivalent to the condition that the eigenvalues 
are real. Unfortunately there seems to be no simple way to then relate this to the 
properties of the function u(t). 
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6 A more general upper bound 



Now let n(t) be an arbitrary everywhere real and nonzero function of t with the 
dimensions of frequency. Then we can rewrite the Schrodinger ODE (pQ) as: 



d 


' Vn ' 




dt 







(tt/tt) 



_ -u 2 (t)/Q ~(n/n) \ I n/Vtt . ' 



n 



(69) 



Again all the matrix elements have been carefully chosen to have the same engineering 
dimension. We can formally solve this in terms of the time-ordered product: 



3? 



{ exp (I 



-cu 2 (i)/n -Un/n) 








B} 





(70) 



to 



The new T matrix is 

t = er 



{ exp u 



-cj 2 (t)/fi -§(fi/fi) 

Note that the matrix T is still real, and that because 
it still follows that T has determinant unity: 



dt 



(71) 







(72) 



a b 
c d 



ad — be = 1. 



(73) 



This means that much of the earlier computations carry through without change. In 
particular as long as at the initial and final times we impose Q(t) — > ujq as t — > tf 
and t — > ti, we still have 



a = ~[a + d + i(b — c)] exp(— icu [tf — U]), 
2 



(3 = - [a — d + i (b + c)] exp(— za;o[£/ + ij]), 

|/3| 2 = ltr{TT r -/}, 

|a| 2 = -tr{TT T + /}. 
4 J 



(74) 

(75) 
(76) 
(77) 
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Now consider the quantity 
X(t) = T(t) T(t) T = sr 



x 



{ exp (/. 



& { exp 



n 



Un/n) n 



di 



t, 



-oo 2 (t)/n -Un/n) 



di 



-i T 



This now satisfies the differential equation 



dX 
~dt 



km) 



n 



with the boundary condition 



X + X 



;(h/tt) -LU 2 (t)/tt 



and 



A brief computation yields 



X(U) = I, 
ti[X] = a 2 + b 2 + c 2 + d 2 . 



dX 
~dt 



n 



+ 



i(n/n) 

-u 2 {i)/Q -|(Q/Q) 

a 2 + b 2 ac + bd 

ac + bd c 2 + d 2 



a 2 + b 2 ac + bd 
ac + bd c 2 + d 2 

■(n/n) -u 2 {t)/Q 



(ft/ft) (a 2 + b 2 ) + 2fi(ac + 6d) fi(c 2 + d 2 ) - (oo 2 /n) (a 2 

-{u 2 /Q){a 2 + b 2 ) + Q{c 2 + d 2 ) -(2u 2 /Q){ac + 6d) - {Q/Q){c 



b 2 ) 

2 i 



Then taking the trace, there is now one extra term 

dtr W / 2 , .2 2 ,2 



df 



fa 2 + b l - & - d l 



2(ac + bd) 



(78) 



(79) 

(80) 
(81) 



f82] 



(83) 



(84) 



Note that if Q(t) — > u then — > and we recover the ODE of the "elementary" 
bound. In this more general setting we now proceed by using the following facts: 

• As previously we note 

\ac+bd\ 2 + \ad-bc\ 2 = a 2 c 2 +2abcd+b 2 d 2 +a 2 d 2 -2abcd+b 2 c 2 = (a 2 +b 2 )(c 2 +d 2 ), 

(85) 
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which implies 



that is 



\ac + bd\ = ^(a 2 + b 2 )(c 2 + d 2 ) - 1, 



2\ac + bd\ = ^/A(a 2 + b 2 )(c 2 + d 2 )-A. 



(86) 
(87) 



Additionally, we use 



+ b 2 - c 2 - d 2 \ = vV + b 2 + c 2 + d?\ 2 - 4(a 2 + b 2 ){c 2 + d 2 ), 



implying 

| a 2 + & 2 - c 2 - rf 2 | 2 + ( 2 |ac + bd\f = \a 2 + b 2 + c 2 + d 2 \ 2 - 4. (89) 

In particular, combining these observations, this means that we can find an angle 9 
(which is in general some complicated real function of a, b, c, d) such that 



2(ac + bd) = vV + b 2 + c 2 + d 2 \ 2 - 4 sin#, 
a 2 + b 2 - c 2 - d 2 = J\a 2 + b 2 + c 2 + d 2 \ 2 -4 cos#, 



whence 

^1 = ^\ a 2 + b 2 + c 2 + d 2 \ 
But for any real 9 we certainly have the inequality 





sin 9 


n 












+ cos 9 




j 




n 









sin# 





+ cos 9 










2 




2 












+ 






n 








n 









implying 



dtrjXj 
dt 



< ^/\a 2 + b 2 + c 2 + d 2 \ 2 -4 



\ 







2 










r w 2 i 






+ 


q - — 




ft 







Therefore 



dtr[X] 
dt 



< v / tr[X] 2 -4 



\ 







2 






"ft" 










+ 


fn--l 




n 







(90) 
(91) 

(92) 



(93) 



(94) 



(95) 



13 



implying 



dtr[X] 



whence 



dcosh~ 1 (tr[A']/2) 
dt 









2 




















+ 

















'ft' 


2 


In-!*" 


2 








+ 

















(96) 



(97) 



so that 



tr[X] = tr[T T T ] < 2 cosh 




ijj- 



dt 



(98) 



Using the general formulae for \a\ 2 and \/3 2 \ in terms of tr{TT T }, and simplifying, 

we see 



< sinh 2 <J - / 1 Jtt 2 + [ft 2 - u 2 ] z dt 



and 



2 L in 



*/ i 



«| 2 < cosh 2 <j - / — ^ n 2 + [n 2 - u 2 Y dt } . 



(99) 
(100) 



This result is completely equivalent to the corresponding result in [TJ; though again 
note that the derivation is largely independent and that it no longer requires one to 
introduce any "gauge fixing" condition, nor need we introduce any WKB-like ansatz. 
The current proof is much more "direct" , and at worst uses simple inequalities and 
straightforward ODE theory. If we work in the space domain instead of the time 
domain and make the translations Q(t) — > <p'{%), oj(t) k(x), we see 



led < cosh 



x f i 

i \<P 



2 — k 2 ] 2 dx 



(101) 



and 



\(3\ 2 < sinh 2 [1 H j i- V / (^) 2 + [(^) 2 -^ 2 ] 2 dx} . (102) 

This is perhaps physically more transparent in terms of the equivalent transmission 
and reflection coefficients 



transmission 



> sech 



x f i 

i \<P 



(103) 
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and 



R < tanh 2 \ - [ * -i-yV') 2 + [(^) 2 -^ 2 ] 2 dx j . (104) 

(For completeness we mention that reference [1] provides a number of consistency 
checks on these bounds by comparing them with known exact results [7].) 



7 The "optimal" choice of Q(t)? 

What is the optimal choice of Q(t) that one can make? Leading to the most stringent 
bound on the Bogoliubov coefficients? The bound we have just derived holds for 
arbitrary Q(t), subject to the two boundary conditions f2(tj) = u = Q(tf) and the 
overall constraint Q(t) ^ 0. Since sinh and cosh are both convex functions, finding 
the most stringent constraint on \(3\ and \a\ is thus a variational calculus problem 
equivalent to minimizing the action 

s= f f -L yj& + [n 2 -u 2 ] 2 dt. (105) 

Jti l"l 

The relevant Euler-Lagrange equations are quite messy, and progress (at least insofar 
as there is any practicable progress) is better made by using an indirect attack. The 
Lagrangian is 

L = — !- sjo? + -uj 2 }\ (106) 
and so the corresponding canonical momentum can be evaluated as 

dQ \Q\ \Jtl 2 + [Q 2 -u 2 } 2 
From the boundary conditions we can deduce 

TTfe) = — = n(t f ). (108) 

The Hamiltonian is now 

n 2 - {n 2 + [n 2 - u 2 } 2 } m 2 -u 2 } 2 

H = 7T Q - L = 1 =J- = [ J =. (109) 

Jtt 2 + [tt 2 - to 2 } 2 |Q| Jtt 2 + [n 2 - uj 2 } 2 
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Unfortunately the Hamiltonian is explicitly time-dependent [via u(t)] and so is not 
conserved. The best we can say is that at the endpoints of the motion 



H(U)=0 = H(t f ). 
By solving for Q as a function of ir and Q we can also write 



and 



n 



H 



VI - vr 2 Q 2 



(110) 



in: 



;ii2) 



Note that Cl at the endpoints is cannot in general be explicitly evaluated in terms of 
the boundary conditions. 

An alternative formulation which slightly simplifies the analysis is to change vari- 
ables by writing 

= oj exp[0(t)], 
where the boundary conditions are now 



e(ti) = o = e(t f ), 

and the action is now rewritten as 



(113) 
(114) 



S 



'9 2 +u; 2 



^10 g — 10 



-i 2 



dt. 



(115) 



Then, in terms of this new variable we have 



L=\ 6 2 + uo 2 



^>1Q g — 10 

^0 



(116) 



with (dimensionless) conjugate momentum 

dL 



71 



86 



;ii7) 



e 2 + ^l 



e 29 _ e -28 



and boundary conditions 



nfc) = 1 =7i(t f ). 



:ii8i 
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The (non-conserved) Hamiltonian is 



H = n 9 — L 



~;20 



9 2 + uu 2 



?26 



4e- 

w 



26 



subject to 



Inverting, we see 



and 



H(t i )=0 = H(t f ). 



7l 



UJ 



7T Z 



e 2e - 



10 



-20 



H 



■vT 



7T 2 LU 



e 2e - 



-20 



;ii9) 



(120) 



:i2D 



:i22i 



This has given us a somewhat simpler variational problem, unfortunately the Euler- 
Lagrange equations are still too messy to provide useful results. 

Overall, we see that while solving the variational problem would indeed result 
in an optimum bound, there is no explicit general formula for such a solution. In 
the tradeoff between optimality and explicitness, we will have to accept the use of 
sub-optimal but explicit bounds. 



8 Sub-optimal but explicit bounds 

From our general bounds 

\P\ 2 < sinh 2 {- jT' ^y^ + ^-c.fdt} , (123) 

and ^ 

\a\ 2 < cosh 2 |1 J ' pV'^ 2 + t fi2 ~ ^ dt ] ' ( 124 ) 
the following special cases are of particular interest: 

Q = ujq: In this case we simply obtain the "elementary" bound considered above. 

Q = uj: This case only makes sense if uj 2 > is always positive. (Otherwise uj and 
hence Vt becomes imaginary in the "classically forbidden" region; the matrix 
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T then becomes complex, and the entire formalism breaks down). Subject to 
this constraint we find 



|/3| 2 < sinh 2 



and 



| or | 2 < cosh 2 

This case was also considered in pp. 







uj 


dtj 


[i 


/" 






hi 


uj 








to 


dtj 
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Jti 


UJ 





(125) 



(126) 



Q = uj e e : This case again only makes sense if uj 2 > is always positive. Subject 
to this constraint we find 



\l3\ 2 < sinh 2 



2 



uj 2 uj 2e 



2-2e 2-2* 
LOq — UJ 



UJ, 



2-2e 




dt 



(127) 



and 



I a I < cosh 



UJ 2 UJ 2 " 



, ,2-2e 2 -2e 
UJq — UJ 



UJ, 



,2-2e 




dt V . 



128) 



This nicely interpolates between the two cases given above, which correspond 
to e = and e = 1 respectively. 



Triangle inequality: Since \J x 2 + y 2 < \x\ + \y\ we see that 



and 



< sinh 



\a\ < cosh 



1 /"*/ 

1 r l s 
24 



dt + 



dt 







n- 


UJ 2 


dt\ 










2 . 


'u 




Q 








UJ 2 






f" 


Q- 




dtj 


2 


Jti 







(129) 



(130) 



These bounds, because they are explicit, are often the most useful quantities to 
calculate. 
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9 The "interaction picture" 

If we split the function u)(t) 2 into an exactly solvable piece ui e (t) 2 and a perturbation 
uj A {t) 2 then we can develop a formal perturbation series for the transfer matrix T, in 
close analogy to the procedures for developing quantum field theoretic perturbation 
theory in the interaction picture. Specifically let us write 

u e {tf = u e (t) 2 + uj A (t) 2 , (131) 

and 

^ = Q(t) T(t) = [Q e (t) + Q A (t)] T{t). (132) 

Now defining 

T(t) = T e (t)T A (t), (133) 
we shall develop a formal solution for T A {t). Consider 

^ = [Qe(t) + Q A (t)] T e (t) T A (t), (134) 

and compare it with 

A m = d J^l T A{t ) + r. W = QJf) Ut) T A{t ) + 7i( t ) (135) 

Therefore 

^-|^ = {TeC*)- 1 Q A (t) T.(t)} T A , (136) 

whence ^ 

T A (t) = ^exp QT {^(t)- 1 Q A (t) T e (t)} dtj . (137) 

For the full transfer matrix T we have 

T{t) = T e {t) x ^exp ^jf {Te(t)- 1 Q A (t) T e (t)} dt) , (138) 

and we have succeeded into splitting it into an exact piece T e (t) plus a distortion due 
to Q A (t). This can now be used as the starting point for a perturbation expansion. 
(The analogy with quantum field theoretic perturbation theory in the interaction 
picture should now be completely clear.) 
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To develop some formal bounds on the Bogoliubov coefficients it is useful to 
suppress (currently) unnecessary phases by defining 



a 



[a + d + i (b - c)\ 



P = - [a - d + i (b + c)] 



(139) 
(140) 



The virtue of these definitions is that for T = T e T& they satisfy a simple composition 
rule which can easily be verified via matrix multiplication. From T = T e Ta we have 



a b 




c d 





a e a a + b e ca a e &a + b e d& 
c e a a + d e ca c e 6a + d e g?a 



Then some simple linear algebra leads to 

P = a e j3 A + (3 e a* A , 
a = a e «a + Pe Pa> 

But then 

\P\ = \P\ 



a e Pa + P e « A 



< 



a e Pa 



Pe &A 



that is 

or the equivalent 
Similarly 



\P\ < \0C e \ \Pa\ + \Pe\ |a A | 



< yJl + \Pef \PA\ + \Pe 



\P\ = \P\ = a e p A + P e a*A 
that is 



> 



at e Pa 



Pe &*a 



(141) 

(142) 
(143) 

(144) 
(145) 

1 + \Pa?. (146) 
= | \a e Pa\ - \Pe «a| I , (147) 



\0L e Pa \ + \Pe «a| 



\P\ > | K\ \p A \ ~ \Pe\ \a A \ 



or the equivalent 



> 



V^ + I&l 2 \Pa\ - \Pe\ yJl + W 



148] 



(149) 



The benefit now is that one has bounded the Bogoliubov coefficient in terms of the 
(assumed known) exact coefficient p e and the contribution from the perturbation Pa- 
Suitably choosing the split between exact and perturbative contributions to u 2 , one 
could in principle obtain arbitrarily accurate bounds. 
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10 Discussion 



In this article we have re-assessed the general bounds on the Bogoliubov coefficients 
developed in [I], providing a new and largely independent derivation of the key results 
that short-circuits much of the technical discussion in [I]. In particular in the current 
article we do not need to "gauge fix", nor do we need to appeal to any WKB-like 
ansatz to get the discussion started. Furthermore we have seen how to extend the 
bounds in [1] in several different ways. 

Considering the fundamental importance of the questions we are asking, it is 
remarkable how little work on this topic can currently be found in the literature. We 
do not feel that the current bounds are the best that can be achieved, and strongly 
suspect that it may be possible to develop yet further extensions both to the current 
formalism, and to the related formalism originally presented in pp. 

Possible extensions might include somehow relaxing the reality constraint on 
n(t) without damaging too much of the current formalism, a better understanding 
of the variational problem defining the "optimal" bound (thus hopefully leading to 
an explicit form thereof), or using several "probe functions" [instead of the single 
function fi(t)} to more closely bound the Bogoliubov coefficients. 



Appendix: Time ordering 

Time-ordered exponentials are a very convenient trick for formally solving certain 
matrix differential equations. Suppose we have a differential equation of the form 

^- = H(t)U(t), (150) 

where U(t) and H{t) are matrices [or more generally linear operators on some vector 
space] and the matrix H(t) is generally not a constant. [So in particular H(ti) need 
not commute with if^)-] In many settings H(t) will be an anti-Hermitian matrix 
in which case U(t) would be unitary — this is not the situation in the current article 
where the matrix H(t) is real and traceless but non- symmetric. 
If Hit) = Hq is a constant then we have the simple solution 

U(t) =exp[H Q t]U(0). (151) 

If H(t) is a constant then we define the formal process of "time ordering" in terms 
of the exact solution U(t) which we know exists because of standard existence and 
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uniqueness theorems. That is 



U(t) = ST \ exp 



H(tf) dt' 



17(0), 



which is equivalent to 



u(t) ir\o). 



(152) 



(153) 



If we take this as our fundamental definition of time ordering then 



dt 



{ 6XP [/ 



H(t') dt' 



H(t)u(t) u-\o) = H(t) sr 



{ 6XP U 



H(t') dt' 



(154) 



But by basic notions of Taylor series expansion 

-t+At 



& { exp 



H(t') dt' 



{/ + H(t) At + 0[(At) 2 ]} 3- \ exp 



exp [H(t) At] & 



{ exp [/ 



#(*') dt' 



H(t') dt' 

+- 0[(At) 2 ]. 
(155) 



Let us now bootstrap this result into a general limit formula for the time ordered 
exponential integral. For simplicity, split the interval (0, t) into n equal segments 
and evaluate H(t) at the points 



I 

n 



then 



^ jexp jf H(t')dt' 



je [0,71-1], (156) 

exp [F(t n _i) At] exp [#(t n _ 2 ) At] . . . (157) 
... exp [if (t x ) At] exp [H(t ) At] + O 



Alternatively 



& jexp 



#(t') dt' 



lim exp [#(*„_! ) At] exp [#(t n _ 2 ) At] . . . (158) 

n— >oo 

. . . exp [H(tx) At] exp [if (t ) At] . 
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This limiting process should remind you of the way the Riemann integral is defined, 
except of course that the H(ti) need not commute with each other so that the order in 
which the matrix exponentials are multiplied together is critically important. This is 
why the product is called "time ordered" . The parameter t can be any real parameter 
- in differential geometry it tends to be a parameter along a curve, sometimes an 
affine parameter, sometimes even arc length, and the product is then sometimes 
referred to as "path ordered", but in general any old parameter would do. 

Note what happens if for some reason the H(tj) do happen to commute with each 
other. Then for instance 



exp [Hfa) At] exp [H(t ) At] -> exp [{H(ti) + H(t )} At] 



(159) 



a result which is not true unless the matrices commute. Continuing in this vein, 
when the matrices do commute we have 



| exp 



Hit 1 ) dt' \ -> lim exp [{if(t n _i) + H(t n _ 2 ) . . .H{t x ) + H(t Q )} At] . 

(160) 

But now the argument of the exponential on the RHS really is the usual Riemann 
integral, so we have 



{ 6XP U 



Hit') dt' 



exp 



H(t') dt' 



(161) 



That is, the time-ordered integral reduces to the ordinary integral whenever the 
matrices H(t) commute with each other. (You could also derive this directly from 
the original differential equation for U(t).) 

In some specific quantum mechanical settings you are more likely to consider the 
slightly different differential equation 



dUjt) 
dt 



-iH(t) U(t), 



(162) 



where H(t) is now the Hamiltonian operator on an appropriate Hilbert space and U 
is the unitary time evolution operator. Then 



U(t) = ST jexp - i jf H{t') dt' | 



U(0), 



(163) 



but note that there is nothing fundamentally new or different here. 
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